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Abstract -The properties of impurities immersed in a large Fermi sea are naturally described in 
terms of dressed quasiparticles: attractive and repulsive polarons, and dressed molecules. Moti- 
vated by recent experiments on narrow Feshbach resonances, we analyze here how the quasiparticle 
properties are affected by a non-zero resonance range. We find two interesting analytic results. 
For large range, the ground state energy close to resonance is shown to become perturbative in 
the inverse range. In the limit of broad resonance instead, we provide a new Tan's relation linking 
the impurity ground state energy E_i to the number of atoms in its dressing cloud AA*'. As a 
corollary, at unitarity one finds AA^ — —Ei/ep, with ep the Fermi energy of the bath. 
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Introduction. — The study of the static and dynamic 
properties of impurities in liquids and crystals provided us 
With a wealth of information on various condensed mat- 
ter systems. As a celebrated example, we may recall how 
Landau showed that the strong interactions experienced 
jby electrons propagating in a crystal could be described 
in very intuitive terms by introducing dressed quasipar- 
ticles named polarons [T]. The investigation of magnetic 
impurities brought to the discovery of the Kondo effect [2] , 
fundamental to understand the electric resistivity of many 
fnaterials at low temperature, while charged impurities in 
Helium liquids [3] played a key role in evidencing vortex 
lattices in the superfluid. 

. Ultracold gases provide an optimal framework for stud- 
ies of many-body quantum effects |3J[5]. The first realiza- 
tion of polarons in an imbalanced homo-nuclear fermionic 
fnixture was reported in an enlightening experiment [5], 
where the energy and residue of the attractive polaron 
Were measured via radio-frequency (rf ) spectroscopy. This 
and later experiments [SlIZ] confirmed that the dressed 
quasiparticles interact very weakly with each other. The 
polaron effective mass was investigated by means of col- 
lective oscillations in ref. [5], and recently the complete 
spectral response of the system was mapped out in ref. [5] . 
There, by exploiting a novel rf spectroscopy scheme, the 
energy of attractive and repulsive polarons was measured, 
and a broad molecular continuum detected. The coherent 
nature of the polarons was probed by driving Rabi oscil- 
lations, providing a precise measurement of their residue. 
The great majority of measurements on strongly- 



interacting fermions employed broad Feshbach resonances, 
but narrow resonances are becoming of increasing rele- 
vance as a new generation of experiments started investi- 
gating hetero- nuclear fermionic mixtures [Timi2j . where 
most known resonances are coincidentally narrow [T5] . 
From a theoretical point of view, at broad resonances the 
relevant physics is a function of a single parameter, the 
scattering length a. The main aim of this paper is the 
study of the richer physics of narrow resonances, where 
quantities depend strongly on a second parameter, the res- 
onance range R* . 

We wish to investigate the low-energy excitation spec- 
trum of a single impurity J, immersed in a large Fermi 
gas of particles f- In presence of a broad resonance, the 
physical scenario is by now well understood [T3H23]- On 
the BCS side, the impurity experiences a weak attractive 
interaction with the Fermi sea, yielding a thinly dressed 
quasiparticle, the attractive polaron. As the attraction 
grows beyond a critical value, it becomes energetically fa- 
vorable for the impurity to bind tightly to an atom in the 
gas, thereby forming a dressed molecule. When the im- 
purity is very light {m^/mi > 7), it becomes possible to 
form even trimers [H]. At positive energies, there exists 
a second branch of polarons, which experience repulsive 
interactions with the gas and are intrinsically metastable, 
due to the presence of lower-lying excitations. 

To be able to investigate narrow resonances, we em- 
ploy in this paper a many-body two-channel model. In 
the strong coupling (broad resonance) limit, our model 
reduces to earlier single-channel approaches, and there it 
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has been vahdated by an ab-initio Monte-Carlo calculation 
[16] . Since we extend earlier treatments towards a weak 
coupling regime (narrow resonance) where the model be- 
comes perturbative in the inverse range and exactly solv- 
able [23] , we expect the results presented here to be accu- 
rate throughout. 

Many-body scattering at narrow resonances. 

The key ingredient in many-body scattering at low ener- 
gies is the T-matrix T(P, lu), describing the collision of an 
t and a ^ atom, with masses m-|- and m^, reduced mass 
rrir, total momentum P and energy w (we have set h — 1). 

At low energies higher partial waves play a negligible 
role, and we may safely restrict ourselves to s-wave scat- 
tering only, parametrized by the scattering length a. The 
usual single-channel T-matrix then reads 

T-liP,^) = mr/{2TTa) - n(P,t^). (1) 

The propagation of the two particles in the medium be- 
tween scattering events is described by the renormalized 
pair propagator n(P, uj) 



n(P,c.) 



dk 

(2^ 



l-/t(k)-/^(P + k) 2m,' 



iO+ - - ^4,p+k 



fc2 



(2) 

with /^(k) = [exp(/3^CTk) + 1] ^ the Fermi function, and 
Ccrk = k'^ / {^rric) — fj-a the kinetic energy measured from the 
chemical potential. The term 2mr/k'^ included in eq. (O 
ensures the convergence of 11 at large momenta, allowing 
us to express the T-matrix in terms of the physical scatter- 
ing length. It is useful to introduce the energy in the center 
of mass reference frame. Ecu = uj — P'^/[2{m^ + mj,)] -f 
Mt + M-l- vacuum limit where both chemical po- 

tentials vanish, one simply finds n(P,w) = — iA;,.(m,/27r), 
where fc, = ^2m,£'cM is the on-shell relative momentum. 
As we see, single channel models containing only a con- 
tact four-fermion interaction fail to describe the physics 
of a narrow resonance, since upon proper renormalization 
of the coupling constant the physics depends on a sin- 
gle parameter, the scattering length a. To make contact 
with experiments, the latter is usually written in the phe- 
nomenological form 



flbg + arcs = abg[l - AB/(S - Bq)]. 



(3) 



Here abg, AB and Bq denote respectively the background 
scattering length, and the resonance width and center. 

A minimal model describing a narrow resonance for- 
mally contains two channels (or liyperfine states), termed 
respectively open and closed. Two particles in the closed 
channel propagate as a molecule with a magnetic moment 
fic differing from the open channel fio by S^l = ^Lc — > 0. 
The two-channel T-matrix may be written as |25II26] 

-1 



r(p,c.) = 



- n(p,w) 



_2'n:a{EcM) 

where we have introduced the energy-dependent length 

AS 



(4) 



a(£^CM) = flbg 1 



B^Bo- Ecu/Sfi 



(5) 



A detailed derivation of this equation is given in the Ap- 
pendix. Comparing with the single-channel T-matrix, we 
see that the inclusion of the second channel simply re- 
sults in a replacement of the scattering length a with the 
energy-dependent quantity a{EcM), or equivalently in a 
magnetic field shift of the resonance center by EcAi/Sfi, 
towards the BCS (BEC) side for positive (negative) Ecm- 
It is interesting to look closely at the low-energy ex- 
pansion of the corresponding on-shell scattering amplitude 
f{kr) = — r(P, aj)mr/27r . In vacuum, we find 



f-^{kr) = a-^ + ikr + R*kl + 0{krY 



(6) 



and higher terms may be neglected at the low-energies 
of interest here. The range parameter R* is related to 
the usual effective range re by the relation R* = —re/2, 
and is given by R* ~ -Rrcsfl ~ Obg/a]^, with R*^^ = 
l/(2mrabgAi?(5/i). Since a — -\-oo as B ^ Bq , one has 
flbgAB > 0, and the range parameter R* is always posi- 
tive. Equation (jG]) reproduces the low-energy physics ex- 
pected from a generic two-channel model with finite range, 
in agreement with refs. [231 [571 [55] . In particular, various 
resonant models such as a square well, a van der Waals po- 
tential [55] , or a confinement- induced resonance [20] show 
the behavior R* cx as a — >■ 0. 

The pole of the truncated vacuum scattering ampli- 
tude eq. ([5]) yields the two-body binding energy Eb = 
— 1/(2™,^*), where the characteristic size of the molecule 
a, = 2i?*/(^l + 4i?*/a - 1) > interpolates between a 
for R* a and V aR* for R* ^ a ^ a^g] moreover, one 
has a* — > Cbg as a — >■ Obg- 

Limiting behaviors. — Given the large number of 
physical parameters in the problem, it is useful to discuss 
here some analytic limits of the theory. 

To start, we consider the negligible background regime 
where jobgl ^ lorosl, see eq. ([3]). This limit is generally 
obtained sufficiently close to resonance, where the scat- 
tering length diverges. In this limit one finds a = Oros, 
R* = and l/a{EcAi) = a-l+2R*mrEcM. To avoid 
an explicit dependence of our results on Obg, we will re- 
strict ourselves in this paper to consider only the negligible 
background regime, and in the following we will drop the 
subscript "res" to avoid cumbersome notationf0. 

We discuss now the condition under which a resonance 
may be termed broad. For this purpose, in eq. ([S]) we 
require |i?*fc^| <C |a~^ -|- ikr\. In ultracold Fermi systems, 
on general grounds one expects bound quasiparticles with 
energy E ^ Eb on the BEC side, and of order —ep on the 
BCS side, while the repulsive polaron will have E ^ +ep. 
Setting fcj. ^ \/2mrE, to be in the broad resonance regime 
one has to satisfy simultaneously the two conditions: 



i?* < a 



and 



F ■ 



(7) 



^Should instead a^g be positive and large, at the two-body level 
the T-matrix admits an extra deep bound state with energy E < 
— l/(2mra^g), as discussed in ref. |26| . yielding a finite lifetime for 
all quasiparticles considered here. 



p-2 



Polarons and dressed molecules near narrow Feshbach resonances 



As a general rule, every s-wave resonance may be con- 
sidered broad at sufBciently low atom densities and suf- 
ficiently close to its center. Given that R* oc g^^, broad 
resonances are obtained whenever there is a strong cou- 
pling between the open and closed channels. 

Quasiparticle equations. — Wc consider a single 
impurity perturbing in a negligible way the surrounding 
Fermi sea. As such, inside the T-matrix we may set /x^ = 
and fi-f = ep, with ep = kp/2m-f and kp = (GTr^n-j-)^/"^ the 
Fermi energy and momentum of the majority atoms. All 
polaron properties arc then given in terms of the impurity 
self-energy, which in the one particle-hole approximation 
(IPHA) and at zero temperature reads [U] 



Ep(p,£;) 



q<kf 



(8) 



where T is the T-matrix. The energies of the two polaron 
branches are given by the two solutions of the equation 



S± =5R[I]p(p,S±+zO+)]. 



(9) 



The polarons residues and effective masses at p = may 
be obtained as Z± = [1 — 9tj3?(Ep)]^^ and 



mi Z± 



(10) 



where the derivatives are taken at the energy of the cor- 
responding quasiparticlcs. The two- and three-body pro- 
cesses leading to the decay of the quasiparticlcs have been 
discussed in detail elsewhere [HHHIIISIISI] ■ 

The energy Eyi of a dressed molecule with momentum 
p is given in the IPHA by the pole of the kernel K of the 
integral equation [TTHH]: 



E 

k' 



K 



k'pq 



(2) 

kk'pq 



-E 



q' 



Akpq> r(p,o) 

.(1) 
kp 



E, 



(1) 
kp 

A'kpq 



E, 



E 



l^k'pq' 



E, 



(1) 



k'q' ^k'p 

r(p,o) 



r(q + p-k,egt-at) 



E, 



(1) 
kp 



(11) 



where 



E, 



(2) 

kk'pq 



Em - ^fct - Cfc't + Cgt ^ C(k+k'-p-q)i 

and iJ^p = Em — ~ C(k-p)4.- Holes and particle 
momenta are restricted to q < kp < k, k' . 

The polaron equation yields the correct result in the 
limit \kpa\ <^ 1, where the attractive/repulsive polaron 
energies become simply Ezf = l-Kajnir. The dressed 
molecule equation instead becomes exact in the deep BEG 
regime, where it reduces to Em — E\, — ep \- liKaADlmj, 
with a AD the atom-dimer scattering length [rTHTOllH^ . 

To gain further insight on the problem, we consider 
analytic (variational) upper bounds for the quasiparti- 
cle energies at p = 0. For the attractive polaron one 
has E_ < Tnin[E']^'^,Ea]. The Thouless pole for the 
attractive polaron E^^^ is determined by the equation 
T-i(fcj?,£'Z'^) = 0. The other bound Eq is obtained from 
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Fig. 1: Quasiparticle energies at resonance and for equal 
masses, versus the resonance width kpR* ■ The dotted lines 
are the Thouless energies E^"^ and El 



cq. by constraining the hole momentum q to 0. Setting 

j/o = —iEo/ep){mr/m^), one finds 



3yo 



—7r{kpa) ^+7rfci?i?*j/o+2+2y^arctan-y/yo. (12) 



On general grounds, one expects that Ep will be closer to 
Eq in the region where the effective mass of the polaron 
is positive, while ET^^ will be a stricter bound when m* < 
and the molecule has lower energy. An upper bound 
for the molecule energy is instead the molecular Thouless 
pole, given by the solution of T-^{0,E^^) = 0. Setting 
y = —{E'^/ep + l)(mr/m^), one finds 







•R*y 



2 -|- 2-\/^arctan-\/^. (13) 



Quasiparticle properties at resonance. We fo- 
cus here on the case = m^, and at [kpa)^^ ~ 0. In 
fig. [T] we show how the energies of the p = polarons 
and dressed molecule are affected by the narrow charac- 
ter of the Feshbach resonance, parametrized by kpR*. In 
the broad resonance limit kpR* ^ 1 the attractive po- 
laron is the ground state of the system with an energy 
= — 0.61ei?, while the dressed molecule experiences 
strong repulsive interactions with the gas, and lies at en- 
ergy Em ~ 0.05ep, more than one Fermi energy above its 
ground state in the absence of interactions. As the res- 
onance becomes narrower, the two quasiparticle branches 
come closer, touch at kpR* = 1.7, and finally the molecule 
becomes the ground state of the mixture. In the asymp- 
totic limit kpR* ^ 1, particle-hole dressing becomes neg- 
ligible and the two energies approach the analytic upper 
bounds given by the respective Thouless poles, E"^^ and 
E'^. In the vicinity of the resonance, from eq. one 
finds 

m-|- 2 — ■K{kpa)~ 



Em 



1 



(14) 



>.r 2 + nkpR* 

Since E"^ > Em > E^ — ep, and i?f, = at resonance, eq. 
(jl4p gives the first terms of a perturbative expansion of the 
ground state energy Em for |fci?a|~^ <C 1 and kpR* ^ 1. 
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Fig. 2: Inverse effective mass mj_/m* (main fig.) and residue 
Z (inset) of the two polarons at resonance for m-f = m^. 

The repulsive polaron is not a good quasiparticle at res- 
onance in the broad resonance limit, as its decay rate is 
too large when compared to its energy p i^^H^ . Nonethe- 
less, a well-defined repulsive polaron appears at unitarity 
for kpR* ^ 1, as in this limit its decay rate becomes 
negligible when compared to its energy 

The inverse effective mass 1/m* and residue Z of the 
polarons are shown in fig. [21 Leaving the broad resonance 
regime, both quantities decrease for the attractive polaron, 
and increase for the repulsive one. In particular, beyond 
the critical width kpR* the effective mass diverges and be- 
comes negative, a usual signature of the ground state tran- 
sition from a polaron to a molecule |16| . Correspondingly, 
increasing spectral weight is continuously transferred from 
the attractive to the repulsive branch. 




Fig. 3: Quasiparticle energies across the resonance for rrif = 
mi. Polarons (molecules) are depicted by thick (thin) lines. 
Continuous blue lines correspond to kpR* = 5, while dashed 
pink ones are for kpR* = 0. The polaron/molecule crossings 
are marked by dots. The dot-dashed blue line is the analytic 
result for the molecule energy with kpR* =5 from eq. (|14p . 
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Quasiparticle properties through the crossover. 

— The features displayed by the quasiparticles at reso- 
nance are actually quite general, and reproduced through- 
out the BEC-BCS crossover and for a large range of mass 
ratios. Our theoretical results for potassium impurities in 
a lithium gas with kpR* ss 1 were presented in ref. [5], 
and agreed with remarkable accuracy with experimental 
observations. Here we focus on the equal masses case, and 
compare in figs. OH] the polaron properties for kpR* = 5 
to known results from the broad case. On the BEG side, 
the attractive polaron and the dressed molecule follow the 
two-body binding energy E^,, and appear less bound than 
in the broad case. On the BCS side, the repulsive in- 
teraction experienced by the dressed molecule is reduced, 
yielding a dressed molecule increasingly bound deep in the 
BCS regime, and a shift of the polaron/molecule crossing 
in the same direction. Remarkably, eq. (fH)) gives a very 
good approximation to -Em close to resonance even for a 
moderately narrow resonance with fci^i?* = 5. The energy- 
dependence of the T-matrix is important for narrow res- 
onances, as noted below eq. and yields sensible shifts 
of all other quasiparticle properties, such as the polarons' 
inverse effective mass and residue shown in fig. |H 



Fig. 4: Inverse effective mass (main figure) and residue (inset) 
of the two polarons for m-j- = mj,. Lines as in fig. [S] 
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Fig. 5: Critical interaction strength of the polaron/molecule 
crossing as a function of the resonance width for various mass 
ratios m-|-/m_|_. Above (below) the line the ground state is a 
molecule (polaron). Inset: energy Ex of the excitations at the 
crossing. The dots mark the interaction strength and energy 
of the crossing as located in the K-Li mixture of ref. |9]. 



p-4 



Polarons and dressed molecules near narrow Feshbach resonances 



The critical interaction strength (kpax)^^ at which the 
polaron and molecule energies are equal is plotted in fig. [5] 
for various mass ratios. For every interaction strength, 
the dressed molecule becomes the ground state of the mix- 
ture at sufhciently large kpR*. For m^/mi = 6/40 and 
kpR* = 0.95, we find —{kpax)"^ ~ 0-6 in good agreement 
with the value obtained in the imbalanccd K-Li mixture of 
ref . [3] (dots in fig. [5]) . The residual small discrepancy may 
partially be ascribed to finite experimental resolution, but 
especially to non-zero temperature effects: and Em 
intersect at a very small angle, and a minute shift of the 
energy levels can yield a noticeable shift of the crossing. 
The inset of fig. (O shows the energy Ex of the two exci- 
tations at the crossing. In the narrow limit kpR* ^ 1, we 



find (kpUx) 



-{kpR*)mr/m^ and E^ oc {kpR*) , as 



one recovers analytically from the corresponding asymp- 
totics Eq ^ E"^ (as in this hmit E. < Eq < < E'^'^). 



Tan's relation for the size of the dressing cloud. 

— The energy density of a Fermi sea containing a few 
impurities {n^ ^ n-j-) is given by [19) 



(3/5)eF"-t + Eini, 



(15) 



where E^ = mm[E^ , Em] is the ground state energy of 
a single dressed impurity. The self-consistent equation 
£'_ = Y,p{E-) at p = gives the polaron energy, and 
a similar one may in principle be written for the dressed 
molecule energy in terms of a self-energy Sm- We restrict 
ourselves here to the broad resonance case, where both 
self-energies are simply functions of the energy uj and of 
the scattering length a. The derivative of eq. ((T5|) with 
respect to v = 1/a yields Tan's contact density C [551 [M] 



C 



de 



Snm, 



dv ^ \ dv 



as ^ di:dEi 
did dv 



n,Z§, (16) 



where E and Z are the ground state (polaron/molecule) 
self-energy and residue, and derivatives are evaluated at 
the quasiparticle energy. The number AA'^ of particles in 
the dressing cloud of the impurity was shown to be |22j 



AiV = - 



E_i 



Z 5E 



dep ep 2aep dv 
Combining the latter two equations, we find 
C 



E^ - ANep. 



(17) 



(18) 



Since the contact remains finite through the crossover, 
at the unitary point of a broad resonance we find 
AN = -Ejep. 

It can be shown that equation (|18|) is actually a limiting 
case of Tan's pressure relation for the impurity problem 
[35] • Indeed, the pressure of an ideal homogeneous Fermi 
gas is Vo = 2ei?n^/5, and its change AV with a change in 
number of particles An^ is AV = V — Vq = 2ep{Anf)/3. 



T 



open 



bg 




E = 



Fig. 6: Diagrammatic structure of tlie two-channel T-matrix. 
Straight and wavy lines indicate respectively propagation in 
the open and closed channels. 

Since Arif = —niAN [55], and Cq = in absence of im- 
purities, eq. (|18p may be rewritten in the form of Tan's 
pressure relation |33j . 



C ^2 

= -P - -e, 

247rmra 3 



(19) 



where e is the energy density of the ideal gas upon addition 
of a few impurities, as given in eq. (jlSp . 



Conclusions. — In this paper we investigated the 
properties of impurities in a Fermi sea via a many-body 
two-channel model, the minimal formalism which captures 
quantitatively the physics of a Feshbach resonance with 
arbitrary width. We calculated the properties of the rel- 
evant quasiparticles, showing how the critical interaction 
strength for the polaron/molecule crossing evolves as a 
function of resonance width and mass ratio. Moreover, the 
ground state energy close to resonance has been proved to 
become perturbative in the small parameter 1/ {kpR*). 

For broad resonances, we have also derived an alterna- 
tive Tan's pressure relation, linking the impurity energy 
E]^ to the number of atoms in its dressing cloud AiV. 

Results presented here are important to understand on- 
going studies on hetero-nuclear ultracold mixtures, spe- 
cially in view of the recent demonstration that narrow 
resonances provide an ideal framework for the realization 
of an elusive state, the strongly-repulsive Fermi gas. We 
hope that our results will stimulate further investigations 
of interacting quantum mixtures. 

Additional remark: shortly after the submission of this 
manuscript, related results appeared in [37] . 

Appendix: T-matrix for narrow resonances. 

We write the two-channel T-matrix as 

T(P, OJ) = Topc„(P, Uj) + Teloscd(P, ^). (20) 

The corresponding Feynman diagrams are sketched in fig. 
[6l and we will discuss them separately below. 

Topcn describes particles propagating only in the 
open channel, undergoing repeated scattering processes 
parametrized by the background scattering length Obg. In 
the ladder approximation, the infinite series of scattering 
events may be resummcd to give 



r-i„(p,a.) = T,-gi-n(p,c.) 



(21) 
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Here Tbg = 2'Kahg/Tnr, and n(P,a;) is the renormalized 
pair propagator given in eq. 

Tcioscd describes instead scattering events that couple 
two particles in the open channel to the closed channel. 
The closed channel molecule experiences an energy shift 
6 fi{B — Bq) with respect to two particles in the open chan- 
nel, and its propagator reads 



D{P, Lj) = [EcM - 5KB - Bo) - E„oi(P, ^)] 



(22) 



We have assumed here that medium effects arising from 
the non-zero population of the closed channel are negligi- 
ble at all times. This assumption is surely justified when 
considering few impurities in a Fermi there may be 

at most as many closed-channel molecules as impurities. 
The two atoms may perform repeated scattering events 
in the open channel before ending up in the closed chan- 
nel, a process overall described by the "dressed vertex" 
F(P,w) =g[l-rbgn(P,w)]-i, with5 = ^TbgABJ/^ |Si 
the renormalized coupling strength between two atoms in 
the open channel and a molecule in the closed channeH. 

The self-energy Emoi describes a molecule coupling to 
two atoms [g], which propagate in the Fermi sea [11] and 
recombine into a molecule through a dressed vertex [V]: 



i],„oi(P,c^) = .gn(P,^)-i/(P,c^), 



(23) 



Only one of the two vertices in cq. ([25)) has been dressed 
to avoid double-counting. The contributions coming from 
the coupling to the closed channel finally yield 

rclo.cd(P,w) = V{P,L0)D{P,L0)V{P,L0). (24) 



Inserting eqs. ([2T|) and ([24| into (j20[) . one finds the com- 
pact expression for the T- matrix given in eq. (|4l). 
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